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Abstract 

In this paper we study a model of s-wave marginal Fermi liquid supercon- 
ductor at two dimensions. Besides the usual Bardeen-Cooper-Schrieffer(BCS) 
point attractive interaction, the fermions in our system also interact via a 
long range repulsive interaction v(q) = go/o 2 - We find that although su- 
perconductivity is not destroyed by the long-ranged repulsion interaction 
(demonstrated by the existence of Meissner effect), nevertheless the BCS 
pairing function Abcs is strongly suppressed resulting in power-law corre- 
lation < Abcs(x)A-bcs{®) >~ l^l -27 - This observation is explained within 
a general picture of marginal Fermi liquid. 
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Since the discovery of high-T c superconductors, there has been enormous interests in the 
study of non- Fermi liquid theories in two dimensions The normal state of high-T c 

superconductors shows rich physical behaviors which seems to deviate from usual Landau 
Fermi liquid theory predictions. For example, in the optimally doped regime, resistivity 
shows a linear temperature dependence, whereas in underdoped regime, the resistivity shows 
a power-law temperature dependence, and with a pseudo-gap in low energy excitation spec- 
trum. Although there exists a lot of theoretical studies on microscopic theories of non-Fermi 
liquid normal states, at present there exists very few systematic investigation into how the 
superconducting state is affected by non-Fermi liquid behaviours. 

In this paper we study a microscopic model of s-wave marginal Fermi liquid super- 
conductor at two dimensions. The Hamiltonian of our model consists of two parts, 
H = H mf i + H BCS , where 

Hmn = E <k)^l(k)Mk) + £ ( la ) 

k,a \g\<gc 

and 

H BCS = -J d 2 x[A(x,t)ij\(x,t)ij\(x,t) + A*(x,t)^(x,t)^(x,t) - ^\A(x,t)\ 2 ] (lb) 

where e(k) = ^- ~ p(%) = J2a Vo-^Vv^) i s the electron density operator, tp a (x)(tpl(x)) 
are electronic annihilation (creation) operators with spin a =|, j, L 2 is the volume of 
the system, q c « k F is a high momentum cutoff, v(q) = go/l 2 , and A(x,t)(A*(x,t)) 
are Hubbard-Stratonovich fields introduced to decouple the (attractive) electron-electron 
interaction U. U is a coupling constant of an unknown mechanism which leads to electron 
pairing. Note that in usual BCS mean-field theory, A(x,t) — > A BC s = —U < ipi(x)ip^(x) >, 
where Abcs is the mean-field BCS gap. The repulsive interaction potential v(q) among 
the electrons may originate from an U(l) gauge field such as in a t-J model or a 

vortex- vortex interaction where each vortex is combined to one electron ||. 

The Hamiltonian H m fi has been carefully studied |5| where it is known that it describes 
a state with Fermi-liquid like particle-hole and collective excitation spectrums, but with 
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nevertheless a vanishing Fermi surface discontinuity z -> | because of orthogonality 
catastrophe effect associated with induced long-ranged density fluctuations when an ad- 
ditional electron is introduced || (marginal Fermi liquid behaviour). Notice that in the 
absence of interaction term v(q), the system would be a regular s-wave superconductor at 
two dimensions, described by a constant pairing function Abcs i n mean-field theory. Here 
we shall examine whether superconductivity is destroyed by the marginal Fermi liquid be- 
haviour, and if it survives, whether and how its behaviour will be modified. Note that even 
though it is believed that in the high-T c cuprates the electron pairing function has a d-wave 
symmetry, interesting conclusions may still be drawn from studying this simple model as we 
shall see in the following. 

To solve this model Hamiltonian we first introduce another Hubbard-Stratonovich field 
4>(q, t) to decouple the repulsive electron-electron interaction. In terms of <fi and A fields, 
the generating functional of the system can be written as, 

\q\<Qc 



2v{q) 



C(t)) 



(2a) 



where T — > oo, 



C(t) = J d 2 x^{x,t)G-\x,t, [(/)})* (x,t), 



(2b) 



and 



G~ (x,t, 



V 2 

id t + (— + ^)a z - A(x, t) - (f>(x, t), 



where a z is a Pauli matrix, ty(x,t) 



4>(x, t) 



</>(£, t), 



^\(x,t) J 



A(f, t) 



0, A(x,t) 
\A*(x,t), j 



(2c) 



and 



Note that the Lagrangian is now quadratic in electron fields 



^ 0, -<j)(x,t) J 
after the HS bosonic fields A(x,t) and <p(x,t) are introduced. 

We shall solve the problem in two steps: (l)for each configuration of 4>(x, t) field, we first 

solve the self-consistent BCS mean-field equations 



G~ l {x,t, [0])G(x,f,t,f, \<p\) = 5{x-tf)5(t-t'), (3) 
A(x,t; [0]) =-U< i[>i(x,t] [<f>])il>i(x,t] [0]) > , 

approximately using an Eikonal expansion |7|,§] or functional bosonization || approach, 
where < A >q is the expectation value of operator A at fixed configuration <p(x,t), and 
(2)with the approximate Green's function G(x,x',t,t'; [4>]) obtained for each configuration 
<p(x,t), we can construct an approximate action S e ff[(f)] in terms of <f){x,t) field only in the 
Eikonal expansion. The ground state expression value < G(x, t) > and other higher-order 
correlation functions can be evaluated as 

^ _ r , !d<pG(x,x\t,t'- [cj)])e- s ^ ... 
<G(x-x ,t-t) >=< G(x, x,t,t;[(f>] >^= JdfaF^fM ' 

etc. Notice that Eikonal expansion can be used as a systematic expansion technique in 
powers of q c jkp [[F]|| for fermions interacting with scalar potential v(q). In particular, the 
marginal Fermi liquid behaviour can be recovered correctly to leading order in Eikonal 
expansion [|7|,|8|,p!0| . The new ingredient we introduce here is to include the BCS mean-field 
theory in a self-consistent way into the formalism. In the following we shall give some brief 
mathematical outlines on how the Eikonal expansion is carried out. First we review the case 
when the BCS term is absent (normal marginal Fermi liquid). 

In this case the matrix Green's function G(x, x', t, [0]) becomes a number function. In 
the Eikonal expansion, we write G as PJlO| 

G(x, x', t, [</>]) = £ go (k, w ) e -*M-fc-+Q £ (-,t;M)] e ^M'-fc-'+Q £ (-'^';M)] ; ( 5 ) 

where g (k,cu) = (uj — e(k))~ 1 , e(k) = k 2 /(2m) — /z, and the function i; [0]) is obtained 
approximately in an Eikonal expansion. To leading order in q c jkp we obtain 

I3l<ffc« w - — 

For the special case of a 1/g 2 potential in two dimension, it turns out that the infra-red 
singularity which leads to marginal- Fermi liquid behaviour comes from plasma oscillations 



with energy uo ~ Vt Q = ygcpo/m, where p is the electron density PJTOl. As a result one 
may neglect the q.k/m term in as far as the leading infra-red singularity is concerned, 
and we obtain 

G(x, x\ t, [<j>]) = G (x -x',t- t 'yiQ{^)-Q{x' ,t')) ^ (g) 

where Q(x, t) ~ Sq> e'^*" 1 ^ ; G (x, t) is the usual free fermion Green's function. 
An Eikonal expansion can also be performed for the effective action S e ff ((/>). To leading order 
in q c /k F , S e ff[(p) is just the RPA effective action 

A1 L q<q c ,oj 

where Xo(q, w) is the usual Linhard function. This result can be understood by noting 
that (i)higher order terms in power of involves sums over extra momentum q < q c , and 
generate higher powers in q c /kp, and (ii)the singular behaviour at small q which is present 
in quadratic term will be absent in higher order terms. Notice also that the marginal 
Fermi liquid behaviour originates from electron-plasmon coupling in the present model pi. 
Therefore, as far as the marginal Fermi liquid behaviour is concerned, the effective action 
may be further simplified by keeping only the plasma-pole contributions [§] to S e f /[</>]. This 
is in contrast to many other models of marginal Fermi liquids where the marginal Fermi 
liquid behaviour is generated from coupling of electrons to anomalous low energy particle- 
hole excitations @||. 

With Eq. (|^), it is easy to show that the corresponding matrix Green's function in the 
presence of BCS attractive interaction which keeps the leading infra-red singularity has form, 



G{x,x',t,t'; 
and with corresponding pairing function 



/ r<BCS i(Q(x,t)-Q(x" ,t')) rpBCS c i(Q(x,t)+Q(x' £)) \ 



T7iBCS* -i(Q(x,t)+Q(3',t')) r< B CS* -i{Q{x,t)-Q{x',t')) 



(?) 



A(f,t;[0]) = ^E^- 2 ^) (8) 
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where G* cs = G BCS (x -x',t- t') and F BCS = F BCS (x -x',t~ t') are the usual BCS 
normal and anomalous Green's functions in the absence of </> field, respectively, with usual 
Fourier transform, 



(iu) 2 - E(k) 

F BCS (kM ~ A 



(ilu) 2 - E(k) 2 ' 

where E(k) = A 2 + e(/c) 2 , A is determined by the usual BCS mean-field equation 1 = 
2^2 J2k ir ( see ©)• We shall assume A << Q in the following discussions. 

As in the normal state, to leading order in q c /kp, the effective action S e ff[4>] of the system 
in the presence of BCS attractive interaction is the quadratic action of superconductors for 
scalar potentials, 



s eff [(t>} = j^ Yl 

<?<|<? c j,n 



2v(q) 



-xo(?,n) + r( g ,n) 



TV 1 \- ^ k v k e q e{-e{k)) 

xo(i, n) = ^ E ( ^ 2 J^^+^ m + ^ - 2 ^ 2 + 4(?))1- 

where e fc (g) = ■ g/m, £ g = q 2 /(2m), u\ = (1 + e(k)/E(k))/2, v\ = (1 - e(k)/E(k))/2, 
and #(a;) is the step function. In the long wave-length limit g — > 0, xo(g ; ^) and T(g, f2) 
are reduced to X o(q -> 0,f2) ~ ^(£) 2 > and r (<7 ^ 0, ft) ~ ^j(^) 4 ~ 0, where p(A) = 
Sfc w fe = Po- F° r a w (<?) = 9d/q 2 potential, the effective action becomes in the limit q — > 
and f2 >> kpqjm^ 



where fto — y9DPo/ m is the plasma frequency. Notice that as far as plasma oscillation is 
concerned, the effective action ^//[^j is unaffected by presence of supeconductivity. This is 
a direct manifestation of Higg's mechanism. 

Using the effective action fllTD , we can calculate the thermodynamical average and cor- 
relation function of the pairing function and the normal and anomalous Green's functions 
by using their expressions in (0) and (|]), respectively. We obtain 



where 7 



, and correspondingly, 



< A(x,t, [0])A*(x',t', 



< A(x,t,[0])A*(fU [0])>^ 



A scs =< A(£, £,[</>]) >^ 




00. 



(13) 



Notice that as the size of the system L 2 — > 00, the pairing function Abcs goes to zero because 
of the long range interaction potential v(q), and its (equal-time) correlation function shows a 
power-law behavior! The normal electron Green's function can also be calculated similarly. 
We obtain 



Notice that a similar expression of the one-electron Green's function is also obtained in the 
normal state with Gq CS replaced by the usual free fermion Green's function. 

There are two effects to be noticed. First of all, we find that the BCS pairing function 
is strongly suppressed by the long-range interaction potential, which seems to indicate that 
superconductivity is suppressed. Secondly, the one-electron Green's function is strongly 
suppressed by a term (g c L)~ 7 at frequency uj < Q , resulting in a 'pseudo-gap' like structure 
in the tunneling density of state for uj < Qq, both in the normal state || and supercon- 
ductivity state. Physically, this 'pseudo-gap' like structure indicates that in the presence 
of long-range repulsive interaction, it is essentially impossible to add (or subtract) an elec- 
tron without creating density fluctuations in the system ||. Notice, however that although 
the differences between superconducting and normal states in the normal and anomalous 
one-electron Green's functions are strongly suppressed by v(q) and are hardly visible when 
A << Qq, nevertheless an order parameter A determined by the BCS mean-field equation 



G(x -x',t- t') ~ G^ cs (x -x / 7 t-t')< 




(14) 



still exists. We shall look into the meaning of nonzero A more carefully by examining directly 
the Meissner effect in the following. 

Following standard linear-response theory analysis, we find in the Coulomb gauge 

J(x, t) = [ d 2 x'P{x -x',t- t')A(x', t') - ^A(x.t) (15) 
J m 

where A(x, t) is the external (transverse) gauge field and J(x, t) is the corresponding induced 
current. The response function P(x,t) is given by 

P(f _ f , t _ t ') = < p(x, x' , t, a- [0]) > 4> , (i6) 

P(x , x', t, t>; [0]) = -^Tr [(V x .V x ,)G(x, x\ t, t'; [<j>])G{? , x, t', t; [0])] , 

where (V.V^G^ = (V x G 1 ).(V x/ G 2 ) + (V x ,G 1 ).(V x G 2 )-(V x .V x ,G 1 )G 2 -G 1 (V x .V x ,G 2 ). 
Using Eq.(7) and the fact that VQ(x, t) introduces extra power of q in Q and is thus higher 
order in g c /^F, we find that to leading order in q c /kF, the Q(x,t) factor does not contribute 
and P(x,t) is given by the usual BCS expression. In particular, the fourier transform 
P(q — > 0, uj = 0) is given by 

Pig :QQ) 1 T E(k)EW-e(k)e(k>)-A* 

which vanishes when A ^ indicating that the system exhibits Meissner effect as in usual 
BCS superconductors. 

The apparent contradiction between existence of Meissner effect and vanishing of BCS 
order parameter Abcs ~ * can be understood in the marginal Fermi liquid picture if we 
assume that the pairing amplitude of quasi-particles A qp = A is nonzero in the ground state. 
The bare fermion pairing amplitude is given by the usual Fermi liquid relation A BC s ~ zA qs 
which is nonzero for usual Fermi liquids with renormalization factor z ^ 0. However, 
z ~ (g c L) -7 and vanishes in the L — > oo limit in the present case and A BC s — > 0. Never- 
theless superconductivity still exists because A qp ^ 0! Notice that this is a rather general 
consequence of a marginal Fermi liquid superconducting state where pairing occurs between 
quasi-particles. In particular, the argument is completely independent of the symmetry of 



the superconducting order parameters. Notice that quantum Monte-Carlo studies on t — J 
model have indicated that the pairing function of electrons Abcs seems to vanish as the 



size of the system increases [nj|. This result was interpreted as indicating the absence of 



superconductivity in t — J model JET]] . Our finding here suggests an alternative possibility: 
superconductivity may still exists in t — J model if the pairing quasi-particles have vanishing 
wavefunction overlap with bare electrons in the system. 

Summarizing, we have studied in this paper the influence of long range repulsive electron- 
electron interaction on an s-wave superconductor at two dimensions. We find that the 
bare electron BCS pairing function Abcs is strongly suppressed by the plasma oscillation, 
with however superconductivity (Meissner effect) kept intact. We interpret the system as 
a marginal Fermi liquid superconductor where the pairing objects are quasi-particles with 
vanishing wavefunction overlap with bare electrons in the system. In particular, we point 
out that this situation might have occurred in the t — J model where the (d-wave) bare 
electron pairing function seems to vanish when size of system increases ||11|| . We also point 
out that a pseudo-gap like structure in both the normal and superconducting states appears 
in the tunneling density of states of bare electrons, as a result of the vanishing wavefunction 
overlap 2 — > 0. An important property of our microscopic model is that the marginal Fermi 
liquid effect comes from coupling of electrons to (high energy) plasmons and is not destroyed 
by superconductivity in contrast with many other models @FJ]. The anomalous behaviour 
we observed in supeconducting state is a direct consequence of survival of marginal Fermi 
liquid effect in the superconducting state and our calculation demonstrates explicitly the 
existence of this new class of marginal Fermi liquid superconductors. 

This work is supported by HKUGC through RGC grant HKUST6143/97P. 
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